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The numerical renormalization group method is used to investigate zero temperature phase tran- 
sitions in quantum impurity systems, in particular in the particle-hole symmetric soft-gap Anderson 
model. The model displays two stable phases whose fixed points can be built up of non-interacting 
single-particle states. In contrast, the quantum phase transitions turn out to be described by inter- 
acting fixed points, and their excitations cannot be described in terms of free particles. We show 
that the structure of the many-body spectrum of these critical fixed points can be understood using 
renormalized perturbation theory close to certain values of the bath exponents which play the role 
of critical dimensions. Contact is made with perturbative renormalization group calculations for 
the soft-gap Anderson and Kondo models. A complete description of the quantum critical many- 
particle spectra is achieved using suitable marginal operators; technically this can be understood as 
epsilon-expansion for full many-body spectra. 
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I. INTRODUCTION 

Zero-temperature phase transitions in quantum impu- 
rity models have recentl y at tracted considerable interest 
(for reviews see Refs. [ll l2l3fl . These transitions can be 
observed in systems where a zero-dimensional boundary 
with internal degrees of freedom (the impurity) interacts 
with an extended bath of fermions or bosons. Examples 
of impurity models with non-trivial phase transitions in- 
clude extensions of the Kondo model where one or two 
magnetic impurities couple to fermionic baths ^, the spin- 
boson model describing a two-level system coupling to 
a dissipative environmenli^, as well as so-called Bose- 
Fermi Kondo models for localized spins interacting with 
both fermionic and bosonic baths. Impurity phase tran- 
sitions are of relevance for impurities in correlated bulk 
systems (e.g. superconductors^) , for multilevel impuri- 
ties like FuUerene molecules^, as well as for nanodevices 
like coupled quantum dots^ or point contacts under the 
influence of dissipative nois o^i^^ . In addition, impurity 
phase transitions have been argued to describe aspects 
of so-called local quantum criticality in correlated lat- 
tice systems. Here, the framework of dynamical mean- 
field theory is employed to map, e.g., the Kondo lattice 
model onto a single-impurity Bose-Fermi Kondo model 
supplemented by self-consistency conditions, for details 
see Ref . [nl 

Diverse techniques have been used to investigate impu- 
rity phase transitions, ranging from static and dynamic 
large- calculationai^, conformal field theory^, pertur- 
bative renormalization group (RG)^'®'^^ and the local- 
moment approach^"*' to various numerical methods. In 
particular, significant progress has been made using the 
numerical renormalization group (NRG) technique, origi- 



nally developed by Wilson for the Kondo problem^^. The 
NRG combines numerically exact diagonalization with 
the idea of the renormalization group, where progres- 
sively smaller energy scales are treated in the course of 
the calculation. NRG calculations are non-perturbative 
and are able to access arbitrarily small energies and tem- 
peratures. Apart from static and dynamic observables, 
the NRG provides information about the many-body ex- 
citation spectrum of the system at every stage of the RG 
flow. Thus, it allows to identify fixed points through their 
fingerprints in the level structure. A detailed understand- 
ing of the NRG levels is usually possible if the fixed point 
can be described by non-interacting bosons or fermions 
- this is the case for most stable fixed points of impurity 
models, e.g., the strong-coupling (screened) fixed point 
of a standard Kondo model. Intermediate-coupling fixed 
points, usually being interacting, have a completely dif- 
ferent NRG level structure, i.e., smaller degeneracies and 
non-equidistant levels. They cannot be cast into a free- 
particle description, with the remarkable exception of the 
two-channel Kondo fixed point which is known to have 
a representation in terms of free Majorana fermionsii. 
In general, the NRG fixed-point spectrum at impurity 
transitions is fully universal, apart from a non-universal 
overall prefactor and discretization effects. 

The purpose of this paper is to demonstrate that a 
complete understanding of the NRG many-body spec- 
trum of critical fixed points is actually possible, by uti- 
lizing renormalized perturbation theory around a non- 
interacting fixed point. In the soft-gap Anderson model, 
this approach can be employed near certain values of 
the bath exponent which can be identified as critical di- 
mensions. Using the knowledge from perturbative RG 
calculations, which yield the relevant coupling constants 
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being parametrically small near the critical dimension, 
we can construct the entire quantum critical many-body 
spectrum from a free-fermion model supplemented by a 
small perturbation. In other words, we shall perform 
epsilon-expansions to determine a complete many-body 
spectrum (instead of certain rcnormalized couplings or 
observables) . Vice versa, our method allows to identify 
relevant degrees of freedom and their marginal couplings 
by carefully analyzing the NRG spectra near critical di- 
mensions of impurity quantum phase transitions. 

The paper is organized as follows. In Sec. ^ we give a 
brief introduction to the physics of the soft-gap Ander- 
son model and its quantum phase transitions. Sec. IIIII 
summarizes the recent results from perturbative RG for 
both the soft-gap Anderson and Kondo models. Section 
llVl describes the numerical renormalization group (NRG) 
approach which is used here to obtain information about 
the structure of the quantum critical points. The main 
part of the paper is Sec. in which we discuss (i) the 
numerical data for the structure of the quantum critical 
points and (ii) the analytical description of these interact- 
ing fixed points close to the upper (lower) critical dimen- 
sion r = (r = 1/2). The main conclusions of the paper 
are summarized in Sec. IVII where we also mention other 
problems for which an analysis of the type presented here 
might be useful. 

II. SOFT-GAP ANDERSON MODEL 

The Hamiltonian of the soft-gap Anderson modeliS is 
given by: 

a 

+ 5]£,cLc,, +VJ2 {PaCka + A J,) • (1) 
ka ka 

This model describes the coupling of electronic degrees 
of freedom at an impurity site (operators f^"^) to a 
fermionic bath (operators c'^l) via a hybridization V. 
The /-electrons are subject to a local Coulomb repul- 
sion U, while the fermionic bath consists of a non- 
interacting conduction band with dispersion Sk- The 
model eq. (Q) has the same form as the single-impurity 
Anderson model^^ but for the soft-gap model we require 
that the hybridization function A(a;) — ttV"^ 5{Ljj—ek) 
has a soft-gap at the Fermi level, A(cl') = A|a;|'", with 
an exponent r > 0. This translates into a local con- 
duction band density of states p(w) = Poller low en- 
ergies. In the numerical calculations we used a band 
where this power law extends over the whole band- 
width D, i.e., from w — ~D/2 to -l-D/2, and we have 
po — (2/1?)''+^ (r 4- l)/2. However, the universal low- 
temperature physics to be discussed in the following does 
not depend on the details of the density of states at high 
energies, and consequently we will use the low-energy 
prefactor of the density of states, pa, to represent the 
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FIG. 1: T = phase diagram for the soft-gap Anderson 
model in the p-h symmetric case (solid line, U — 10~^, £/ = 
—0.5 • 10~^, conduction band cutoff at -1 and 1) and the p-h 
asymmetric case (dashed line, £/ = —0.4 • 10"'^); A measures 
the hybridization strength, A{uj) = A\ui\^ . 

dimensionful energy scale of the problem. Assuming a 
particle-hole symmetric band, the model is particle- 
hole symmetric for ej = —U/2. 

The soft-gap Anderson model with < r < oo 
displays a very rich behaviour, in particular a continuous 
transition between a local-moment (LM) and a strong- 
couphng (SC) phase. Figure ^ shows a typical phase 
diagram for the soft-gap Anderson model. In the particle- 
hole (p-h) symmetric case (solid line) the critical coupling 
Ac di verges at r = ^, and no screening occurs for r > ^ 
(Refs. Il9l20j) . No divergence occurs for p-h asymmetry 
(dashed line)i^. 

We now briefly describe the properties of the fixed 
points in the soft-gap Anderson and Kondo modelsiS. 
Due to the power-law conduction band density of states, 
already the stable LM and SC fixed points show non- 
trivial behaviourii2i2&. The LM phase has the properties 
of a free spin i with residual entropy Simp — fc_Bln2 
and low-temperature impurity susceptibility Ximp — 
l/{4kBT), but the leading corrections show r-dependent 
power laws. The p-h symmetric SC fixed point has very 
unusual properties, namely S'imp = 2rfcBln2, Ximp = 
r/ (SksT) for < r < i. In contrast, the p-h asymmetric 
SC fixed point simply displays a completely screened mo- 
ment, S'imp — Tximp — 0- The impurity spectral function 
follows a oj^ power law at both the LM and the asym- 
metric SC fixed point, whereas it diverges as w"*" at the 
symmetric SC fixed point - this "peak" can be viewed 
as a generalization of the Kondo resonance in the stan- 
dard case (r = 0), and scaling of this peak is observed 
upon approaching the SC-LM phase boundarjii^iSS. At 
the critical point, non-trivial behaviour corresponding to 
a fractional moment can be observed: S'imp = kBCs{r), 
Ximp = C^{r) / [ksT) with 0$, being universal func- 
tions of r (see Refs. IT^l2ll) . The spectral functions at 
the quantum critical points display an u;"*" power law 
(for r < 1) with a remarkable "pinning" of the critical 
exponent. 
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III. RESULTS FROM PERTURBATIVE RG 

The Anderson model is equivalent to a Kondo 
model when charge fluctuations on the impurity site 
are negligible. The Hamiltonian for the soft-gap Kondo 
model can be written as 



H 



(2) 



where s(0) = Z^fcfc'crcr' 4<T'^crcr'Cfc'cr'/2 is the conduction 
electron spin at the impurity site r = 0, and the con- 
duction electron density of states follows a power law 
p{ijj) = pol'^r ^ above. 



A. RG near r = 

For small values of the density of states exponent r, 
the phase transition in the pseudogap Kondo model can 
be accessed from the weak-coupling limit, using a gener- 
alization of Anderson's poor man's scaling. Power count- 
ing about the local moment fixed point (LM) shows that 
dim [J] = — r, i.e., the Kondo coupling is marginal for 
r = 0. We introduce a renormalized dimensionless Kondo 
coupling j according to 



PoJ = fJ. ''j 



(3) 



where /i plays the role of a UV cutoff. The flow of the 
renormalized Kondo coupling j is given by the beta func- 
tion 



(4) 



For r > there is a stable fixed point at j* = cor- 
responding to the local- moment phase (LM). An unsta- 
ble fixed point, controlling the transition to the strong- 
coupling phase, exists at 



J 



(5) 



and the critical properties can be determined in a double 
expansion in r and j^. P-h asymmetry is irrelevant, i.e., 
a potential scattering term E scales to zero according to 
/?(e) = re (where pqE = p^^e), thus the above expansion 
captures the p-h symmetric critical fixed point (SCR). As 
the dynamical exponent l/v = r + 0{r'^), diverges as 
r — !■ 0^ , r = plays the role of a lower-critical dimension 
of the transition under consideration. 



B. 



RG near r = 1/2 



For r near 1/2 the p-h symmetric critical fixed point 
moves to strong Kondo coupling, and the language 
of the p-h symmetric Anderson model becomes more 
appropriate^^. First, the conduction electrons are inte- 
grated out exactly, yielding a self-energy = V^Gco for 



the / electrons, where Gco is the bare conduction elec- 
tron Green's function at the impurity location. In the 
low-energy limit the / electron propagator is then given 
by 



Gf{iujn) ^ = ii^n - iAo sgn(w„) \ujnY' 



(6) 



where the \uJnY self-energy term dominates for r < 1 
and the prefactor Aq is 



(7) 



Equation © describes the physics of a non-interacting 
resonant level model with a soft-gap density of states. 
Interestingly, the impurity spin is not fully screened for 
r > 0, and the residual entropy is 2rln2. This precisely 
corresponds to the symmetric strong-coupling (SC) phase 
of the soft-gap Anderson and Kondo modeliS,. 

Dimensional analysis, using dim[/] = (1 — r)/2 [where 
/ represents the dressed fermion according to eq. 
now shows that the interaction term U of the Anderson 
model scales as dim[[/] — 2r — 1, i.e., it is marginal at 
r = 1/2. This suggests a perturbative expansion in U 
around the SC fixed point. We introduce a dimensionless 
renormalized on-site interaction u via 



(8) 



The beta funcion receives the lowest non-trivial contri- 
bution at two-loop order and readai^ 



I3{u) = (1 - 2r)u 



3(7r-21n4) 



u'^ + Oiu^). (9) 



For r < 1/2 a non-interacting stable fixed point is at u* — 
- this is the symmetric strong-coupling fixed point, it 
becomes unstable for r > 1/2. Additionally, for r < 1/2 
there is a pair of critical fixed points (SCR, SCR') located 
at u*^ = 7r2(l - 2r)/[3(7r-21n4)], i.e.. 



±4.22 VI72" 



(10) 



These fixed points describe the transition between an un- 
screened (spin or charge) moment phase and the symmet- 
ric strong-coupling phasei^. 

Summarizing, both Q and ^ capture the same criti- 
cal SCR fixed point. This fixed point can be accessed ei- 
ther by an expansion around the weak-coupling LM fixed 
point, i.e., around the decoupled impurity limit, valid for 
r ^ 1, or by an expansion around the strong-coupling SC 
fixed point, i.e., around a non-interacting resonant-level 
(or Anderson) impurity, and this expansion is valid for 
1/2-r < 1. 



IV. NUMERICAL RENORMALIZATION 
GROUP 

Here we describe the numerical renormalization group 
method, suitably extended to handle non-constant cou- 
plings A(w) (see Refs. 11912011221ft . This method aUows a 
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non-perturbative calculation of the many-particle spec- 
trum and physical properties in the whole parameter 
regime of the model eq. in particular in the low- 
temperature limit, so that the structure of the quantum 
critical points is accessible, as discussed in Sec. 

A detailed discussion of how the NRG can be applied 
to the soft-gap Anderson model can be found in Refs.ll9l 
Here we focus on those aspects of the approach 
necessary to understand how information on the fixed 
points can be extracted. 

The NRG is based on a logarithmic discretization of 
the energy axis, i.e. one introduces a parameter A and 
divides the energy axis into intervals [— A~", — A~("+-'^^] 
and [A-("+i), A'"] for n = 0,1,... .,00 (see Refs. lll2l . 
With some further manipulations the original model can 
be mapped onto a semi-infinite chain with the Hamilto- 



(T 



oc 



+ lcr + C„^i(jCn(j 



with 



/ dwA(cj) 



(11) 



(12) 



For a p-h symmetric conduction band, all the on-site en- 
ergies Sn vanish. If, in addition, the power law in A(w) 
extends up to a hard cut-off tOc (we set lUc — 1), an ex- 
act expression for the hopping matrix elements t„ can be 



givei 



tn = A-"/2 



1 l-A-('^+2) 



r + 2 1 - A-('^+i) 
^^ -1/2 r 

X 1 _ A-(2"+-+i) 1 



1 - A"("+'^+^) 



A-(2«+r+3)l '/Ji3) 



for even n and 

t„ = A-("+'-)/2 

I _ ^-(2n+r+l 



r -f 1 1 - A-('^+2) 

r + 2 1 - A-('-+i) 
-1/2 r 



1 - A"("+^) 



I _ ^-(2n+r+3) 



-1/2 



(14) 



for odd n. The semi-infinite chain is solved iteratively by 
starting from the impurity and successively adding chain 
sites. As the coupling i„ between two adjacent sites n 
and n+1 decreases as A^"/^ for large n, the low-energy 
states of the chain with n + 1 sites are generally deter- 
mined by a comparatively small number Ng of states close 
to the ground state of the n-site system. In practice one 
retains only these iVg states from the n-site chain to set 
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FIG. 2: Flow diagrams for the low-energy many-body ex- 
citations obtained from the numerical renormalization group 
for the three different fixed points of the p-h symmetric soft- 
gap Anderson model (exponent r — 0.4). N is the number 
of iterations of the NRG procedure, A the NRG discretiza- 
tion parameter. Solid lines: {Q,S) — (1,0), dashed lines: 
(Q,5) = (0,1/2). 

up the Hilbert space for the n+1 site chain, thus prevent- 
ing the usual exponential growth of the Hilbert space as 
n increases. Eventually, after tt-nrg sites have been in- 
cluded in the calculation, addition of another site will not 
change significantly the spectrum of many-particle exci- 
tations; the spectrum is very close to that of the fixed 
point, and the calculation may be terminated. 

In this way, the NRG iteration gives the many-particle 
energies En for a sequence of Hamiltonians Hn which 
correspond to the Hamiltonian eq. (|ll|l by the replace- 
ment 



N-l 



(15) 



An example for the dependence of the lowest lying en- 
ergy levels on the chain length (the flow diagram) is 
given in Fig. ^ for the soft-gap Anderson model with 
r = 0.4, £1 = 2, U/D = 10~^ and A = 0.0075; the 
parameters used for the NRG calculations are A = 2 
and Ns = 300. The states are labelled by the quantum 
numbers Q (which characterizes the number of particles 
measured relative to half- filling) , and the total spin, S 
[solid hues in Fig. [2 are for {Q,S) — (1,0), dashed lines 
for {Q,S) — (0,1/2)]. As mentioned above, the energy 
scale is reduced in each step by a factor A^^^. To allow 
for a direct comparison of the energies for different chain 
lengths, it is thus convenient to plot A^^^E^ instead of 
the eigenvalues En of the A^-site chain directly. Note 
that here and in the following we use the convention that 
the energies shown in the flow diagrams are proportional 
to the bandwidth D. 

As is apparent from Fig. the properties of the sys- 
tem in this case do not change further for chain lengths 
JT-NRG > 120. Without going into details now, one can 
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FIG. 3: Dependence of the many-particle spectra for the three 
fixed points of the p-h symmetric soft-gap Anderson model on 
the exponent r: SC (black dot-dashed lines), LM (blue dashed 
lines), and the (symmetric) quantum critical point (red solid 
lines) . The data are shown for the subspace Q = 1 and 5 = 
only. 



show that the distribution of energy levels for > 120 
in Fig. [2t is characteristic of the SC phase of the model 
(see SecEJ. 

If by contrast we choose instead a value of A = 0.006, 
we obtain the flow diagram shown in Fig. Here it 
is evident that the fixed point level structure is entirely 
different from the SC solution, and indeed this particu- 
lar pattern is now characteristic of the LM phase of the 
model. We can thus conclude, simply from inspection 
of the two flow diagrams, that the critical Ac separating 
the SC and LM phases of the soft-gap Anderson model 
for the model parameters specified, lies in the interval 
[0.006,0.0075]. 

Tuning the value of A to the critical value Ac, results 
in the flow diagram of Fig.|2h- Apparently, the structure 
of the fixed point at Ac neither coincides with the SC 
nor with the LM fixed point. It is clear that it cannot be 
simply constructed from single-particle states as for the 
SC and LM fixed points. An important observation is 
that certain degeneracies present in the SC or LM fixed 
points are lifted at the QCP. As shown in the following 
section, a further hint on the structure of the QCPs is 
given by the dependence of their many-particle spectra 
on the bath exponent r. 



V. STRUCTURE OF THE QUANTUM 
CRITICAL POINTS 

In Fig. O the many-particle spectra of the three fixed 
points (SC: dot-dashed lines, LM: dashed lines, and QCP: 
solid lines) of the symmetric soft-gap model are plotted 
as functions of the exponent r (for a similar figure, see 
Fig. 13 in Ref . ITol) . The data are shown for an odd num- 
ber of sites only and we select the lowest lying energy 
levels for the subspace Q = 1 and S = 0. 



FIG. 4: Difference AE between the energy levels of QCP and 
SC fixed points close to r = 1/2 in a double-logarithmic plot. 
The inset shows the values of the exponents obtained from a 
fit to the data points. 

As usual, the fixed point structure of the strong cou- 
pling and local moment phases can be easily constructed 
from the single-particle states of a free conduction elec- 
tron chain. This is discussed in more detail later. Let us 
now turn to the line of quantum critical points. What 
information can be extracted from Fig.|31to understand 
the structure of these fixed points? 

First we observe that the levels of the quantum critical 
points, Eff^QCpii"), approach the levels of the LM (SC) 
fixed points in the limit r ^ (r 1/2): 

lim {£:Ar.QCp(r)} = {EN,LMir ^ 0)} , 
lim {EN,QCp{r)} = {^w,sc(r= 1/2)} , (16) 

where {. . .} denotes the whole set of many-particle states. 

For r —> 0, each individual many-particle level 
En, Qcpir) deviates linearly from the levels of the LM 
fixed point, while the deviation from the SC levels is pro- 
portional to -\/l/2 — r for r 1/2. This is illustrated 
in Fig. 01 where we plot a selection of energy differences 
AE between levels of QCP and SC fixed points close to 
r = 1/2. The inset shows the values of the exponents 
obtained from a fit to the data points. For some levels, 
there are significant deviations from the exponent 1/2. 
This is because the correct exponent is only obtained in 
the limit r — > 1/2 (the QCP levels have been obtained 
only up to r = 0.4985). 

Note that the behaviour of the many-particle levels 
close to r = 1/2 has direct consequences for physical 
properties at the QCP; the critical exponent of the local 
susceptibility at the QCP, for example, shows a square- 
root dependence on (1/2 — r) close to r = 1/2, see Ref.lTgl 

In both limits, r — > and r 1/2, we observe that de- 
generacies due to the combination of single-particle lev- 
els, present at the LM and SC fixed points, are lifted at 
the quantum critical fixed points as soon as one is mov- 
ing away from r = and r = 1/2, respectively. This 
already suggests that the quantum critical point is inter- 
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acting and cannot be constructed from non-interacting 
single-particle states. 

In the remainder of the paper we want to show how to 
connect this information from NRG to the perturbative 
RG of Sec. mil We know that the critical fixed point 
can be accessed via two different epsilon-expansionsSii^ 
near the two critical dimensions, and, combined with 
renormalized perturbation theory, these expansions can 
be used to evaluate various observables near criticality. 
Here, we will employ this concept to perform renormal- 
ized perturbation theory for the entire many-body spec- 
trum at the critical fixed point. To do so, we will start 
from the many-body spectrum of one of the trivial fixed 
points, i.e., LM near r = and SC near r = 1/2, and eval- 
uate corrections to it in lowest-order perturbation theory. 
This will be done within the NRG concept working di- 
rectly with the discrete many-body spectra correspond- 
ing to a finite NRG chain (which is diagonalized numer- 
ically). As the relevant energy scale of the spectra de- 
creases as A~"/^ along the NRG iteration, the strength 
of the perturbation has to be scaled as well, as the goal 
is to capture a fixed point of the NRG method. This scal- 
ing of the perturbation follows precisely from its scaling 
dimension - the perturbation marginal at the value of r 
corresponding to the critical dimension. With the proper 
scaling, the operator which we use to capture the differ- 
ence between the free-fermion and critical fixed points 
becomes exactly marginal [see eqs. H21II and (|37|l below]. 



As usual, the structure of the fixed point spectra depends 
on whether the total number of sites is even or odd. To 
simplify the discussion in the following, we only consider 
a total odd number of sites (the flow diagrams of Fig. [21 
arc all calculated for this case). For the LM fixed point, 
this means that the number of sites, -I- 1, of the free 
conduction electron chain is even, so A^ in eq. (|17|l is odd. 

The single-particle spectrum of the free chain with an 
even number of sites, corresponding to the diagonalized 
Hamiltonian 

Hc,N = X! ''P Ccr^P<^ J (18) 

is sketched in Fig. El (Note that the gp have to be eval- 
uated numerically.) 



El P=1 
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A. Perturbation theory close to r = 

Let us now describe in detail the analysis of the devia- 
tion of the QCP levels from the LM levels close to r = 
(the case r = 1/2 is discussed in Sec. IV B]) . An effec- 
tive description of the LM fixed point is given by a finite 
chain with the impurity decoupled from the conduction 
electron part, see Fig. jSl 

12 N 




o-o-o- - -o 

free conduction electron chain 



decoupled impurity 



FIG. 5: The spectrum of the LM fixed point is described 
by the impurity decoupled from the free conduction electron 
chain. 

The conduction electron part of the effective Hamilto- 
nian is given by 

N-l 

HcM = ^" {4icrCn+la + cl+i^C„^) . (17) 

an—O 



FIG. 6: Single-particle spectrum of the free conduction elec- 
tron chain eq. 11811 . The ground state is given by all the levels 
with p < filled. 

As we assume p-h symmetry, the positions of the 
single-particle levels are symmetric with respect to with 

ep = -6_p , p= l,3,...,Af , (19) 

and 

E ■ (20) 

p P— — N, p odd 

Note that an equally spaced spectrum of single-particle 
levels is only recovered in the limit A ^ 1 (see Fig. 6 in 
Ref.lH for the case r — 0. 

The RG analysis of Sec. Illll tells us that the critical 
fixed point is perturbative accessible from the LM one 
using a Kondo-type coupling as perturbation. We thus 
focus on the operator 

K = air)fiN) Sin,p ■ so , (21) 

with the goal to calculate the many-body spectrum of 
the critical fixed point via perturbation theory in H'j^ for 
small r. The function a{r) contains the fixed-point value 
of the Kondo-type coupling, and f{N) will be chosen 
such that H'js^ is exactly marginal, i.e., the effect of H'j^ 
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on the many-particle energies decreases as K~^/'^ which 
is the same N dependence which governs the scahng of 
the many-particle spectrum itself. The scaling analysis 
of Sec, mil eqs. ^ and ||SJ), suggests a parametrization 
of the coupling as 



a{r) 



Po 



ar . 



(22) 



where po is the prefactor in the density of states, and /i 
is a scale of order of the bandwidth - such a factor is re- 
quired here to make a a dimensionless parameter. Thus, 
the strength of the perturbation increases linearly with 
r at small r (where fj.~^/po = D + 0{r) for a featureless 
density of states). 
The qualitative influence of the operator 5'imp ■ sb on 
the many-particle states has been discussed in general in 
Ref. for finite r and in Refs. for r — 0. Whereas 

an antiferromagnetic exchange coupling is marginally rel- 
evant in the gapless case (r = 0), it turns out to be 
irrelevant for finite r, see Ref. This is of course con- 
sistent with the scaling analysis of Sec. IIIII the operator 
(|21|l simply represents a Kondo coupling, with a tree-level 
scaling dimension of dim[J] = — r. A detailed analysis of 
the A^-dependence of the operator 5'inip ■ so shows that it 
decreases as A~^^^'^A~^/'^ = A~^('"+^^/^ with increasing 
A^. Consequently, we have to choose 

f{N) ^ A^'-/2 . (23) 

This result also directly follows from dim[J] — ~r: As 
the NRG discretization yields a decrease of the running 
energy scale of A~^/^, the S'imp-s'o term in lf?T)) scales 
as A~^^/'^. The function /(A) is now simply chosen to 
compensate this effect; using eq. (1^ the operator H'^ 
becomes exactly marginal. 

Now we turn to a discussion of the many-body spec- 
trum. The relevant ground state of the effective model 
for the LM fixed point consists of the filled impurity level 
(with one electron with either spin "f or |) and all the con- 
duction electron states with p < filled with both | and 
I, as shown in Fig. |^ Let us now focus on excitations 
with energy ei -f 63 measured with respect to the ground 
state. Figure [7| shows one such excitation; in this case, 
one electron with spin [ is removed from the p = —3 
level and one electron with spin [ is added to the p = 1 
level. The impurity level is assumed to be filled with an 
electron with spin "f , so the resulting state has Q = and 
Sz = -1-1/2. In total, there are 32 states with excitation 
energy ei +63. These states can be classified using the 
quantum numbers Q, S, and Sz- 

Here we consider only the states with quantum num- 
bers Q = 0, 5 = 1/2, and Sz = 1/2 (with excitation 
energy ei -I- e3) which form a four-dimensional subspace. 
As the state shown in Fig.[7|is not an eigenstate of the to- 
tal spin S, we have to form proper linear combinations to 
obtain a basis for this subspace; this basis can be written 
in the form: 



1^1) 



V2 



H 



-3 



FIG. 7: One possible excitation with energy ei -I- £3 and quan- 
tum numbers Q = Q and Sz = +1/2. 



IV'2) = 



Jg/|(eiV-3T-eU-3i) + ^/|eU-3x 



IV'3) = -i=/|(et^e-iT+4?-u) iV^i 



IV^4) = 



1^0) 



(24) 



where the state |?/'o) is given by the product of the ground 
state of the conduction electron chain and the empty im- 
purity level: 



IV'o) 



.P<0 



® |0)imp 



(25) 



The fourfold degeneracy of the subspace (Q = 0, 
5"= 1/2,5*2 = 1/2) of the LM fixed point at energy £1-1-63 
is partially split for finite r in the spectrum of the quan- 
tum critical fixed point. Let us now calculate the influ- 
ence of the perturbation H'^ on the states . . • |^/'4), 
concentrating on the splitting of the energy levels up to 
first order. Degenerate perturbation theory requires the 
calculation of the matrix 



{ipi\H'j^\i).j) , «,j = i. 



(26) 



and a subsequent calculation of the eigenvalues of {Wij} 
gives the level splitting. 

Details of the calculation of the matrix elements Wij 
are given in Appendix ^ The result is 



{W,,} = a(r)/(A) 






V3 
4 ' 








4 ' 



















^7 








4 ' 


-i/3 



(27) 



with 7 = [|aoiP - |ao-3n and /? = [|aoiP + |ao-3p]- 
The A^-dependence of the coefficients aop (which relate 
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r=0.03 



the operators coa and ^p^, see eq. ^A8\i ) is given by 

|aop|'ocA-^'-/2A-^/2 , (28) 
(see also Sec. Ill A in Ref-El). Numerically we find that 



7 = -0.1478 •A-^'-/2a-^/2 
f3 = 2.0249 • a-^'-/2A-^/2 , 



(29) 



where the prefactors depend on the exponent r and the 
discretization parameter A (the quoted values are for r = 
0.01 and A = 2.0). The matrix {Wyj^^Qg^ then takes 
the form 



y ir=0.01 



= a(r = 0.01)A-^/2 

-0.064 
-0.064 -1.013 

-0.064 











-0.064 -1.013 



(30) 



Diagonalization of this matrix gives the first-order cor- 
rections to the energy levels 

A£;i(r = 0.01) = AS3(r = 0.01) 

= a(r = 0.01) A-^/2 . 0615) 
A£;2(r = 0.01) = A£^4(7- = 0.01) 

= a(r = 0.01) A"^/2 . 0.0004 (31) 

with 

EN,QCp{r = Om,i) = EN,LMir = 0.01,i) 

+ AE,{r^0.01) , (32) 

(j = 1, . . . 4). Apparently, the fourfold degeneracy of the 
subspace {Q ~ 0, S ~ 1/2, 5*^ = 1/2) with energy ei -I- 63 
is split in two levels which are both twofold degenerate. 

We repeated this analysis for a couple of other sub- 
spaces and a list of the resulting matrices {VFy } and the 
energy shifts AE is given in Appendix IXI 

Let us now proceed with the comparison of the pertur- 
bative results with the structure of the quantum critical 
fixed point calculated from the NRG. For our specific 
choice of the conduction band density of states, the rela- 
tion (|22|) yields a{r) = ar D for small r (where fi^ ~ 1). 
Using the corresponding equations for the energy shifts in 
Appendix ^ we observe that a single parameter a must 
be sufficient to describe the level shifts in all subspaces, 
provided that the exponent r is small enough so that the 
perturbative calculations are still valid. A numerical fit 
gives a « 1.03 for A = 2.0, (the A-dependence of a is 
discussed later, see Fig.|5)|. 

Figure |S| summarizes the NRG results together with 
the perturbative analysis for exponents r close to 0. A 
fiow diagram of the lowest lying energy levels is shown 
in Fig. for a small value of the exponent, r = 0.03, 
so that the levels of the QCP only slightly deviate from 



a) 



< 

z 
LLI 




QCP 
LM 



b) 



< 

z 
LU 



— QCP 
-- LM 
+ pert. 



50 

N 



0.03 



0.07 



FIG. 8: a) Flow diagram of the lowest lying energy levels for 
r — 0.03; dashed lines: flow to the LM fixed point; solid lines: 
flow to the quantum critical fixed point, b) The deviation of 
the QCP levels from the LM levels increases linearly with r. 
This deviation together with the splitting of the energy levels 
can be explained by the perturbative calculation (crosses) as 
described in the text. 



those of the LM fixed point. As discussed above, the 
deviation of the QCP levels from the LM levels increases 
linearly with r, see Fig. |5Jd. We indeed get a very good 
agreement between the perturbative result (crosses) and 
the NRG-data (lines) for exponents up to r « 0.07. The 
data shown here are for the subspaces {Q — 0, S = 1/2, 
= 1/2) and energy 2ei (the levels at E^k^^'^ « 1, 
see Appendix m|l and (Q = 0, 5 = 1/2, 5*^ = 1/2) and 
energy ei -l-ea (the levels at SjvA^^^ ~ 2, see the example 
discussed in this section). 

In the NRG, the continuum limit corresponds to the 
limit A ^ 1, but due to the drastically increasing nu- 
merical effort upon reducing A, results for the contin- 
uum limit have to be obtained via extrapolation of NRG 
data for A in, for example, the range 1.5 < A < 3.0. 
Figure El shows the numerical results from the NRG cal- 
culation together with a linear fit to the data: a (A) — 
0.985 -t- 0.045(A - 1.0). Taking into account the increas- 
ing error bars for smaller values of A, the extrapolated 
value a{K ^ 1) « 0.985 is in excellent agreement with 
the result from the perturbative RG calculation, which 
is directly for the continuum limit and gives a = 1.0. 



B. Perturbation theory close to r = 1/2 

To describe the deviation of the QCP levels from the 
SC levels close to r = 1/2, we have to start from an effec- 
tive description of the SC fixed point. This is given by a 
finite chain including the impurity site with the Coulomb 
repulsion J7 = at the impurity site and a hybridiza- 
tion V between impurity and the first conduction elec- 
tron site, see Fig. ^1 

Note that the SC fixed point can also be described by 



9 



1.5 



1.0 



0.5 



0.0 



1.0 



^-e— e-=e- e- 



1.5 



2.0 

A 



2.5 



3.0 



respect to with 

eo = , £/ = 
and 



Z = 2,4,...,(iV + l) 



l=N+l 



E 



E 

/=-(Af+l), I oven 



(35) 



(36) 



The ground state of the effective model for the SC fixed 
point is fourfold degenerate, with all levels with I < 
filled and the level I = either empty, singly (f or I) or 
doubly occupied. 



FIG. 9: Dependence of the coupling parameter a on the NRG- 
discretization parameter A. The circles correspond to the 
NRG data and the solid line is a linear fit to the data: q(A) — 
0.985 + 0.045(A- 1.0). 



the limit V —I- oo and finite U which means that impurity 
and first conduction electron site are effectively removed 
from the chain. This reduces the number of sites of the 
chain by two and leads to exactly the same level structure 
as including the impurity with U — 0. However, the 
description with the impurity included (and U = 0) is 
more suitable for the following analysis. 



U=0 



1 



N 



■o-o-o- 



o 



FIG. 10: The spectrum of the SC fixed point is described by 
the non-interacting impurity coupled to the free conduction 
electron chain. 

The corresponding effective Hamiltonian is that of a 
soft-gap Anderson model on a finite chain with N + 2 
sites and e f = U — (i.e., a p-h symmetric resonant 
level model). 



sc,JV 



= V 



[n 



COo 



c,N 



(33) 



with Hc^N as in eq. ((TTIl . 

As for the effective description of the LM fixed point, 
the effective Hamiltonian is that of a free chain. Fo- 
cussing, as above, on odd values of TV, the total number 
of sites of this chain, N + 2, is odd. The single-particle 
spectrum of the free chain with an odd number of sites, 
corresponding to the diagonalized Hamiltonian 



E 



(34) 



is sketched in Fig. ^2 As we assume p-h symmetry, the 
positions of the single-particle levels are symmetric with 



1=4 
1=2 



£^=0 



1=0 



H 



l=-2 
l=-4 



FIG. 11: Single-particle spectrum of the free conduction elec- 
tron chain eq. 13411 . The ground state is fourfold degenerate 
with all the levels with I < filled and the level I = either 
empty, singly (| or J.) or doubly occupied. 

According to Sec. IIIII the proper perturbation to ac- 
cess the critical fixed point from the SC one is an on-site 
repulsion, thus we choose 



(37) 



(n/o- = flfa) with the strength of the perturbation 
parametrized as 



(3{r)=pi'^-'plv'f3^/l/2^ 



(38) 



see SecHTTI Note that pl{r = 1/2) = 9/(2^3) for a fea- 
tureless power-law density of states with bandwidth D. 
The N dependence of the operator — i) (n/| — i) 
is given by a('-i/2)^A-^/2 = A^'-'i)^, so we have to 
choose 



fiN) = A(i/2--)w 



(39) 



This again follows from the scaling analysis of Sec. IIIII 
the on-site repulsion has scaling dimension dim[f/] = 2r— 
1. Thus the (n/| — i) — i) term in H'j^ (|?7|l scales 
as A-^^'^~^/'^\ and f{N) H39|l compensates this behavior 
to make H'j^ exactly marginal. 
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We continue with analyzing the low-lying many-body 
levels. Similar as above, we focus on one specific example, 
these are excitations with energy 2e2 measured with re- 
spect to the ground state and quantum numbers Q — —1, 
S — 0, and Sz — 0. This subspace is two-dimensional and 
the basis is given by 



IV-i) 

1^-2) 



^fe-2T+ek-2i) IV'o) , (40) 



with 



(41) 



(Note that in this definition of |V'o), the I = 0-level is 
empty.) 

The two-fold degeneracy of this subspace is lifted for 
r < 1/2 in the spectrum of the quantum critical points. 
The matrix Wij = {t/jilH'j^l^jjj) = 1,2) is given by 



72| 



2-2k + k2 
2V2k 



(42) 

with K = |a/oP/|a/2p. The A^-dependence of the coef- 
ficients \afi\ (which relate the operators and ^i^, see 
eq. (|B4|) ') is given by 

la/H' oc A('■-lW^ (43) 

Numerically we find that 

\af2\^ = 0.0366 • (i?/F)2A('^-i)^/2 

|a/oP = 0.0930 • (i?/F)2A('^-i)^/2 , (44) 

where the prefactors depend on the exponent r and the 
quoted value is for r = 0.499. The matrix {Wij}^^i^^gg 
then takes the form 



499 



I3{r = QA99){D/V)^A~^/^ 

0.0044 0.0094 
0.0094 0.011 



(45) 



Diagonalization of this matrix gives the first-order cor- 
rections to the energy levels 



{i = 1,2). We repeated this analysis for a couple of other 
subspaces and a list of the resulting matrices {Wij} and 
the energy shifts AE is given in Appendix IbI 

The comparison of the pcrturbative results with the 
numerical results from the NRG calculation is shown in 
Fig. 112b . As for the case r « we observe that a single 
parameter /3 is sufficient to describe the level shifts in 
all subspaces, provided the exponent r is close enough 
to r — 1/2 so that the pcrturbative calculations are 
valid. For A = 2.0 we find /3 « 70 and the A ^ 1 
extrapolation results in /3(A — > 1) w 73.0 ± 5.0 (the 
error bars here are significantly larger as for the extra- 
polation of the coupling a) . The results from pcrturba- 
tive RG, Sec. mil specifically eqs. ^ and ifTUI) . yield 
P{r) = /^2'--Vo^^27r2w*. This gives (3 = 83.3. 

Similar to Fig. |S1 above, we show in Fig. 112b a flow 
diagram for an exponent very close to 1/2, r = 0.4985, 
so that the levels of the QCP only slightly deviate from 
those of the SC levels. As discussed above, this devia- 
tion is proportional to \/l/2 — r, see Fig. 112b. The data 
shown here are all for subspaces with {Q = —1, S = 0, 
Sz = 0); the unperturbed energies E of these subspaces 
are: 



AEi{r 
AE2{r 

with 



0.499) 
0.499) 



P{r = 0.499) (i:'/V')*A-^/2 • (-0.0023) 
f3{r = 0.499) (i?/y)4A-^/2 • (0.018) 

(46) 



En, QCP {r 
EN,sc{r 



0.499, i) 
0.499,?) 



E = 0: the levels at SjvA^/^ ^ ggg App.IBT 
E = £2: the levels at EnA^^^ « 0.8, see App.|B 



• E = 2e2- the levels at E^A^/"^ w 1.6, see the ex- 
ample discussed in this section, 



AEi{r = 0.499) , (47) 



• E = ei: the levels at EnA^/"^ k. 1.8, see App.IbH 

• E = 362: the levels at EnA^/"^ w 2.4. 

We again find a very good agreement between the pcr- 
turbative results (crosses) and the NRG data (lines). 

Thus we can summarize that our renormalized pertur- 
bation theory for the NRG many-body spectrum works 
well in the vicinity of both r = and r — 1/2. In princi- 
ple, from the many-body spectrum (and suitable matrix 
elements) all other observables like thermodynamic data 
and dynamic correlation functions can be calculated. We 
note that the convergence radius of the epsilon-expansion 
for the levels seems to be smaller than that of the di- 
rect epsilon-expansion for certain observables like criti- 
cal exponents and impurity susceptibility and entropy, 
see Ref . Ill 



VI. CONCLUSIONS 

Using the quantum phase transitions in the soft-gap 
Anderson model as an example, we have demonstrated 
that epsilon-expansion techniques can be used to deter- 
mine complete many-body spectra at quantum critical 
points. To this end, we have connected information from 
standard pcrturbative RG, which yields information on 
critical dimensions and parametrically small couplings, 
and from NRG for the many-body spectra of free-fermion 
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r=0.4985 



a) 





- SC 

QCP 
+ pert. 
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0.495 
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APPENDIX A: DETAILS OF THE 
PERTURBATIVE ANALYSIS AROUND THE 
LOCAL MOMENT FIXED POINT 



12: a) Flow diagram of the lowest lying energy levels 
r — 0.4985; dashed lines: flow to the SC fixed point; 



FIG 

for 

solid lines: flow to the quantum critical fixed point, b) The 
deviation of the QCP levels from the SC levels is proportional 
to •y/l/2 — r. This deviation together with the splitting of the 
energy levels can be explained by the perturbative calculation 
(crosses) as described in the text. 



fixed points. Together, these can be used to perform 
renormalizcd perturbation theory for many-body spec- 
tra of interacting intermediate-coupling fixed points. For 
the soft-gap Anderson model, which features two lower- 
critical dimensions at r = and r = 1/2, correspondingly 
two different approaches can be utilized to capture the 
same critical fixed point: Near r = a Kondo term has 
to be added to a free-fermion chain with a decoupled im- 
purity, whereas near r = 1/2 an on-site repulsion is used 
as a perturbation to the non-interacting Anderson (or 
resonant-level) model. These perturbations lift the large 
degeneracies present in the non-interacting spectra, and 
accurately reproduce the critical spectra determined in 
NRG calculations at criticality. 

Vice versa, our method will be useful in situations 
where the effective low-energy theory for the critical point 
is not known: a careful analysis of the many-body spec- 
trum near critical dimensions yields information about 
the scaling dimension and structure of the relevant oper- 
ators. 

For instance, a plot similar to Fig.|21can be calculated 
for the spin-boson model, using the numerical renormal- 
ization group method as in Ref. 0. Preliminary results 
(not shown here) indicate that the many-particle levels 
of the QCP approach the levels of the delocalized (local- 
ized) fixed point in the limit s — )■ {s 1), with s the 
exponent of the bath spectral function J(w) cx cj*. 

We envision applications of our ideas to more complex 
impurity models, e.g., with two orbitals or two coupled 
spins, as well as to non-equilibrium situations treated 
using NRGS4. 



In this Appendix, we want to give more details for the 
derivation of the matrix Wij eq. (|27|l which determines 
the splitting of the fourfold degeneracy of the subspace 
{Q ^ 0, S = 1/2, 5*2 = 1/2) of the LM fixed point at 
energy ei -I- £3. We focus on the matrix element Wi2'. 

W12 = {i^i\H'jy\^P2) = a{r)f{N) (V'l IS'imp • so|^2) . (Al) 

The strategy for the calculations can be extended to the 
other matrix elements and the other subspaces, for which 
we add the perturbative results at the end of this ap- 
pendix without derivation. The operator S'imp • sq is de- 
composed in four parts: 



5*1, 



■ So 



1 

V 

2 ^ im 



Q+ t 



1 



r'S'i, 



1 



1 



(c^ 



2 -'imp \^'^0T''''T 

SO that W12 can be written as 

W12 = a{r)f{N) ^ [I + n + III - IV] 



with 



I 



(A2) 



(A3) 



(A4) 



and the other terms accordingly. With the definitions of 
and IV'2) of eq. (|24|l we have: 



I = 



-3i 



Ho) . 

(A5) 



With 5;^p = /|/| we immediately see that the terms 
containing /-f'5'imp/| drop out. The remaining impurity 
operators, f^S^^fJ^, give unity when acting on \tpo) so 



one arrives at 



x/3 



[la + lb] 



(A6) 



12 



with 

la = (V'olClai^TcJxCoT^i^'f-sil'/'o) 

lb = (V-oielsiCucJ^coTeiV-sil^o) . (A7) 

To analyze la and lb, the operators c^]. have to be ex- 
pressed in terms of the operators : 

COa ^^aop'£.p'a , 4<T = X! "Op^'J ' (A8) 

p' p 

with the sums over p and p' defined in eq. (j2()|l . This 
gives 

la = ^a5paop'(V'o|Cl3T?iT^pi^p'TCnC-3ilV'o) ■ (A9) 
pp' 



Following similar arguments as above one obtains 



The only non-zero matrix elements of eq. (|A9|) are for 
p = p' = -3: 

la aS-3"O-3(^o|Cl3|ClT^-3i^-3T^ItC-3il'0o) 

- -|«o-3|' . (AlO) 
Similarly, the term lb gives 

lb = ^aSpaop'(?Ao|Cl3x6iCpxCp'TClT^-3il'/'o) 
pp' 

= laoil' , (AH) 
so that in total: 



1= -= h|ao-3|' + |aoiP] 



V3 



(A12) 



III 



The next term II = {'^i\S\^pc\^-^CQi\tp2) gives zero due to 
the combination of impurity operators: fifff^--- with 
/l from (t/jiI and /|/| = 5i7„p. 

The third term III = (V'l |5f^pcJ^cot |-02) gives 

X (eiV-3T - ei\e-3i) l^o) , (A13) 

where the term with --2=y|^t^^_3| from \1lj2) has already 
been dropped. So we are left with four terms 

III = [Ilia - Illb + IIIc - Illd] , (A14) 

with 

Ilia = (?Ao|d3|6T/T'S'impcS|Cot/|^{^C-3TlV'o) , 

Illb = (V'o|^l3|6T/T'5'impcJ|Cot/|CuC-3ilV'o) , 

IIIc = (V'o|d3i6i/T'S'topcJxCoT/|.ei\C-3TlV'o) , 

Illd - (V^0|d3i6i/T5impc5TC0T/|ei\C-3ilV'0) ■ 



Ilia 



aopi 



(A16) 



where the p in takes the values 

p= 1,-1, -5, -7,... -TV 



then 



and 



Illb = IIIc = 



Illd = ij^'VopI' 



where the p in X]p" takes the values 

p= -1,-3,-5,-7,...- 
This gives for the third term 
1 



(A17) 



(A18) 



III 



[Ilia - Illd] 



1 



4\/3 
1 

4Vl 



L P 



[jaoi 



|ao-3|' 



(A19) 



The calculation of IV proceeds very similarly to III and 
one obtains 



III = -IV 
so that we finally arrive at 



12 = «(r)/(A^)i(|aoi|2-|ao_3p) 



(A20) 



1 1 

^ + + 2^ 
\/3 4V3 



- a(r)/(A^)iV3(|aoi|2-|ao-3|') 



(A21) 



(A15) 



We performed a similar analysis for a couple of other 
subspaces. Here we list the results from the perturba- 
tive analysis for three more subspaces together with the 
corresponding basis states. 



1. Q = 0, S = 1/2, = 1/2, E = 2ei 

This subspace has the same quantum numbers Q, S 
and ^2 as the one discussed above, so that the details 
of the calculation are very similar. The differences origi- 
nate from the position of particles and holes in the single- 
particle spectrum of Fig. This reduces the dimension- 
ality of the subspace from four to two. 
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The corresponding basis can be written as 



m ■ 

(A22) 

The first-order corrections are given by the 2x2 matrix 



{W,,} = a{r)f{N) 



^7 

#7 -¥ 



(A23) 



with 7 = |aoip — jao-iP and /3 = jaoiP + |ao-iP- Due 
to the particle-hole symmetry of the conduction band we 
have |aoi| = |ao-i|; therefore, the off-diagonal matrix el- 
ements vanish and the effect of the perturbation is simply 
a negative energy-shift only for the state \tp2)- 



{W.,} = a{r)f{N) 







|aoi| 



(A24) 



This effect can be seen in the energy splitting of the first 
two low-lying excitations in Fig. |S1 



2. Q = -1, S = 0, E ^ -e_i 

There is only one configuration for this combination of 
quantum numbers and excitation energy: 



IV')-^(/|C-iT+/|^u)l^o) 



(A25) 



APPENDIX B: DETAILS OF THE 
PERTURBATIVE ANALYSIS AROUND THE 
STRONG COUPLING FIXED POINT 

The main difference in the calculation of the matrix 
elements {Wij} for this case is due to the structure of 
the perturbation, see eq. (|37|l . Furthermore the ground 
state of the SC fixed point is fourfold degenerate and the 
perturbation partially splits this degeneracy, as discussed 
in the following. 



1. Q = 0, S=l/2, S, = l/2, E = Q 

This is one of the four degenerate ground states at the 
SC fixed point: 



(Bl) 



with j-^o) defined in eq. H41|) . 

The perturbative correction is given by 

(MKly^i) = lmfiN)il - la/oh (B2) 
which corresponds to the energy shift of the ground state: 
AE,^^P{r)f{N){l-\afo\') . (B3) 

The coefficients aji are defined via the relation between 
the operators f^^ and ^[p: 



^ -1, S = 0, E = 



The first-order perturbation keeps the state in this one- 
dimensional subspace and the energy correction is given 
by 



AE = {ib\H'N\^) - --air)f{N)\ao-i\^ . (A26) 



3. Q = -l, 5 = 0, £ = -e_3 

The difference to the previous case is the position of 
the hole in the single-particle spectrum. The state is now 
given by 



1^) = ^(/|e_3T+/|e-3i)|V'0) , 



V2 

with the energy correction 



(A27) 



AE = im'NW = --a(r)/(A^)|ao-3p ■ (A28) 



This state is also a ground state in the U = case: 



IV'2) = I^O) 



(B5) 



The calculation of the first-order correction for \tlj2) gives 
(HKm = i/3(r)/(iV)(l + |a/oh • (B6) 

This means that the ground state including the effect of 
the perturbation is given by in eq. IjBll) and the state 
IV'2) appears as an excited state. For a comparison with 
the energy levels shown in the NRG flow diagrams, where 
the ground state energy is set to zero in each iteration, 
we subtract the perturbative correction of the ground 
state (AEi ) from the energies of all other excited states. 
Subtracting this energy shift from eq. (|B6|I gives the net 
energy correction for the \ip2) state: 



AE2^Pir)f{N)\afo\' 



(B7) 
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3. = -1, S = 0, _E = £2 
The state corresponding to this subspace is given by: 

iV's) = ^(^^T^-2T + e^xe-2i)ivo) . (B8) 



The first-order correction reads 

'1 



(1- |a/or) + 3|a/or|a/„2|^ 
(B9)^ 

Subtracting the energy correction for the ground state 
results in 

AEs = 3/3(r)/(iV)|a/ona/-2|' ■ (BIO) 

4. = -1, S = 0, _E = £4 
Similarly for the state 

1^4) = ^(f0T^-4T + ^^i^-4i)|V'o) , (Bll) 



the first-order correction is given by 



{4'a\h'm = mm) 



-(l-|a/oh + 3|a/ona/-4|' 
(B12)' 



and subtracting the energy correction for the ground 
state results in: 



A^4 = 3/3(r)/(iV)|a/ona/- 



(B13) 
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